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We have applied a microscopic model for single photon emission in neutral current interactions 
on nucleons and nuclei to determine the number and distributions of such events at the Super- 
Kamiokande detector, for the flux and beam exposure of the T2K experiment in neutrino mode. 
These reactions represent an irreducible background in electron-(anti)neutrino appearance measure¬ 
ments aimed at a precise measurement of mixing angle ^13 and the CP violating phase. We have 
obtained a total number of photon events that is twice larger than the one from the NEUT event 
generator (version 5.1.4.2) used in the analysis of T2K data. Detailed comparisons of energy and 
angular distributions for the and fluxes have also been performed. 


I. INTRODUCTION 


More than 15 years of dedicated experimental studies have established the oscillations of three flavors of massive 
neutrinos. Recently, a new piece in the puzzle has been added by the determination of the so-called reactor neutrino 
mixing angle ^ 13 . The first indication with a 2.5 cr significance of a nonzero value of ^13 was provided by the T2K 
experiment in a study of Ve appearance in a beam |l|. Afterwards it has been precisely measured from Df, 
disappearance in nuclear reactor neutrino experiments [2 t 4- The significance of the T2K z/g appearance result has 
now reached 7.3 a The increasing precision in these experiments creates a window of opportunity to determine the 
CP violating phase in the lepton sector. Indeed, the tension between reactor data and T2K favors a 5cp = — 7 r /2 @ 
at 90% C.L., although the picture is still far from clear because the MINOS combined disappearance and z/g 
appearance prefers a(5cp = 7 r /2 0 . 

Further progress in this direction requires a better control over systematic errors and, in particular, of irreducible 
backgrounds. For this purpose, a better understanding of neutrino interactions with matter is mandatory. The 
ongoing effort in this direction encompasses more precise measurements of different (anti)neutrino cross sections on 
nuclear targets, theoretical work aimed at a better description of weak reactions on both nucleons and nuclei, and 
improvement of the Monte Carlo simulation codes; see the reviews of Refs. for different aspects of these problems. 

Super-Kamiokande (SK), the far detector of the T2K experiment, is a water Cherenkov detector and, as such, is 
incapable of discriminating the diffuse rings of originated in charged current interactions by electron neutrinos 
from those created by photons. The largest part of such a background originates in 7 r° production in neutral current 
(NC) interactions (NCtt*^) when the two photons from —>■ 77 produce overlapping rings or when one of them is not 
observed. Another relevant source is the NC single photon emission (NC 7 ) reaction. Although NCtt^ has a larger cross 
section than NC 7 , the 7 r° background can be reduced with dedicated reconstruction algorithms, while the NC 7 one 
remains irreducible. Indeed, in the latest T2K analysis the NCtt® background was reduced by 69% with respect 
to the previous appearance selection [l^. In this context, the relative relevance of the NC 7 channel is significantly 
enhanced. 

The interest in a detailed theoretical study of the NC 7 reaction EMI followed the observation of an excess of 
electron-like events at low reconstructed energies in the MiniBooNE detector, in both neutrino and antineutrino 
modes SIS- It was suggested that an anomalous contribution to NC photon emission could be responsible for 
this E3' However, in spite of the fact that first studies indicated that NC 7 indeed accounted for the excess E^i more 
recent analyses, considering nuclear effects and realistic acceptance corrections [13, [HI, obtain a number of photon- 
induced electron-like events which is consistent with the estimate made by MiniBooNE, using a poor resonance 
production model tuned to the experiment’s own NCtt^ measurement. It has also been proposed that additional 
photons from electromagnetic heavy neutrino decays could be at the heart of the MiniBooNE anomaly [13, [13| , which 
would have implications for other experiments such as T2K and MicroBooNE. 

Here, we apply the microscopic model of Ref. [l3|, briefly described in Sec. im to predict the number of NC 7 events 
at the SK detector, as well as their energy and angular distributions, for the flux and beam exposure of the latest T2K 
z/e appearance study [3] . The event number and distributions are compared to those of the NEUT event generator (23| , 
from which the T2K estimates are obtained. 
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II. OUTLINE OF THE THEORETICAL MODEL 

The main features of our model for NC photon emission induced by (anti)neutrinos on nucleons and nuclei are 
presented in the following. A detailed account is given in Ref. [l5l| . with comparisons to previous results. A concise 
overview of the different models can be found in Ref. 


A. NC7 on 

The differential cross section for the reactions 

viik) + N{p) viik') 
vi{k) + N{p) Di{k') 

on nucleons (protons or neutrons) in the laboratory frame is 

d^cF(u,v) _ Ej 
dE^dn{kj) ~ \k\ 167r2 

where k = {E^, k), k' = {E', k') and k^ = {E^, k^)] G denotes the Fermi constant. The leptonic tensor 

+ k',k^ + < 7 ^,^ ± le^ao^pk'^kd , ( 2 ) 

with q = k — k' and the fully antisymmetric tensor defined such that eoi 23 = +1, is contracted with the hadronic one 

= I ^^s\p' + k^-q-p){N^\rmN){N^\fmNr, (s) 

where p = (mjv,0) and p' = {E'j^^p'). A sum (average) is performed over the final (initial) spin states. The matrix 
element of the hadronic current 


nucleons 

+ N{p') + 7(A;^) , 

+ E[{p') + 7(A;^) 

given by 

( 1 ) 

|U| 


(iV7|j^(0)|iV)=n(p')r)(ru(p)e;(fc^), 


(4) 


with e{k~f) being the photon polarization four vector, is 


defined [l5l | by the set Feynman diagrams shown in Fig. [T] 



FIG. 1. Feynman diagrams for the hadronic current of NC photon emission off nucleons considered in Ref. [l^ . The first two 
diagrams represent direct and crossed baryon-pole terms with nucleons and nucleon resonances in the intermediate state: BP 
and GBP with B = N, A(1232), A(1440), A(1520), A(1535). The third diagram denotes t-channel pion exchange: tiEx. 

At neutrino energies around and below 1 GeV, where most of the T2K flux is concentrated (see Fig. [3]), the NC 7 
processes are dominated by the weak excitation of the A(1232) resonance followed by its radiative decay. Non-resonant 
nucleon-pole contributions that close to threshold, are fully constrained by the chiral symmetry are also important. 
We have also taken into account mechanisms with nucleon excitations from the second resonance region [Pn Af(1440), 
Di 3 Ai(1520) and Ai(1535)] as intermediate states. Among them, the Ai(1520) contributes most significantly (see 
Fig 5. of Ref. jH). The pion-exchange mechanism arises from the anomalous vertex and gives a very small 
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contribution to the cross section [T^. Other t-channel p and w exchange amplitudes, not considered in Ref. (TB| . 
coming from the corresponding anomalous vertices [l^, have turned out to be small [ll| , particularly according to 
the latest estimate [^ . 

In terms of the amputated amplitudes introduced in Eq. o 

^ “ 9 -p)Tr [(/ + mN)Tf^P{^ + mNh°{T'^p)'’7°] , (5) 

with 

a = BP, CBP, ttEx, and B = N, A(1232), iV(1440), 7V(1520), A(1535). ( 6 ) 

a 

Explicit expressions for these amplitudes can be found in Ref. [T5j |. They are functions of phenomenological nucleon 
and nucleon-to-resonance vector and axial form factors (partially) constrained by (quasi)elastic and pion production 
reactions induced by electrons and neutrinos on hydrogen and deuterium targets. 


B. NC photon emission in nuclei 

On nuclear targets, NC photon emission can be incoherent, 

vi{k) + Az —>■ ) + 7 (^ 7 ) T ^) 

vi{k) + Az —>■ Di(k') + 7(^7) T ^ 5 

with the hnal nucleus in any excited state or broken, or coherent 

Vl{k) + Az\gs{PA) Vl{k') + Az\gs{p'A) + ) 

Vl{k) + Az\gs{PA) vi{k’) + Az\gs{p'A) + t 

when the nucleus is left in the ground state. Equation m remains valid in both cases, with the hadronic tensor 
replaced by the pertinent one, Wj^coh O'" ^coh- 



FIG. 2. Diagrammatic representation of the one-particle-one-hole-photon (Iplhy) contributions to the self-energy in nuclear 
matter. The black dots represent the Z^N 7 N elementary amplitudes of Fig [T] To obtain the imaginary part, the 
intermediate states intersected by the dashed line have to be placed on the mass shell. 

In the incoherent case. Ref. [l5l| applies a many-body scheme for neutrino propagation in nuclear matter adapted to 
(semi)inclusive reactions on finite nuclei using the local density approximation |27| . In this framework, the hadronic 
tensor is obtained as the imaginary part of the contributions to the Z selfenergy with a single photon in the intermediate 
state. To the lowest order in a density expansion, corresponding to a one-particle-one-hole nuclear excitation (Figure[2]) 
this tensor reads 


Coh 


<p‘;h7 = 0(g°)0('?° 


( 2 mw)^ 


f kj,k^,k^)A’'j^ , 

N=p,n'' ^ 


(7) 


where 


^7 = ^Tr 


-b mAf) 7° + + Wiv) 


( 8 ) 
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in terms of the amputated ZN —>■ 77 V elementary amplitudes introduced in the previous section. This tensor is 
evaluated at an average nucleon hole four momentum p = (p) to simplify the calculation. Our choice for (p^) and 
discussions in support of this approximation are given in Refs. |l5l. l28l| . The definition and explicit expressions for 
the Lindhard function U ^{q — k^,kp ,kp) can be found in Ref. [27|. It depends on the local density of protons or 
neutrons in the nucleus through the local Fermi momenta kp{r) = [37r^p7v(r)]^/^. 

For the coherent reaction, the nucleon wave functions are not altered by the interaction; the amplitudes should be 
summed over all nucleons in the target. This leads to [T^ 


with 


^c'oh = 


SjE^ - g°) 


(A%) 


* 


(9) 


and 


APPiq,kj)= I d\ [p^{t)VPp + Pr.{r)vpp] , 


2Tr[(j^ + m^)7°r 


/ip] 

N \ „0 


-qp 


iik-i-q) 


( 10 ) 


( 11 ) 


in terms of the elementary ZN -p N'y amputated amplitudes of Eq. ([6]) and Fig[TJ In these amplitudes, energy 
conservation is accomplished by imposing q^ = which is justified by the large nucleus mass, while the transferred 
momentum is assumed to be equally shared between the initial and final nucleons. This approximation is discussed 

in Refs. i[il. 

In both incoherent and coherent processes, the broadening of the A(1232) resonance 


Fa —t Fa — 2ImEA(p) (12) 

in the nuclear medium is included. The resonance decay width is reduced to Fa because the hnal nucleon in A —>■ -kN 
can be Pauli blocked but, on the other hand, is increased by many body processes A N ^ N N, AN -p N N tt and 
ANN -p N N N. For these decay channels, accounted in ImEA, we have taken the parametrizations of Ref. [s^- 
This collisional broadening results in a reduction of the NC 7 cross sections in nuclei [1^ . 


C. Error budget 

The theoretical model described above has uncertainties in the treatment of both hadronic interactions and nuclear 
effects. We have performed an error analysis by propagating the uncertainties listed in Table HI assuming that they 
are uncorrelated and Gaussian distributed. At the hadronic level, we consider errors in the leading A —A(1232) axial 
coupling (7^(0) and in the parameter controlling i ts dependence (axial mass M^a) according to the analysis of weak 
pion production on deuterium performed in Ref. |3l| . The uncertainties in the N — A largest helicity amplitudes A 1/2 
and A 3/2 at = 0 are also included by taking the relative errors from the PDG estimates [s^. The unitary isobar 
model MAID [s^ , from which our nucleon-to-resonance electromagnetic form factors are determined [ 1 ^ 0 , does not 
provide these errors. Uncertainties in the mechanisms with N* intermediate states can be large, particularly in the 
axial transition currents which are poorly known. They can nevertheless be neglected because at the low energies of 
the peak in the T2K neutrino flux (Figure [S]), the contribution from these terms is quite small. Finally, in the case of 
the nucleon form factors present in the NP and CNP amplitudes, we neglect errors in the vector form factors and 
axial coupling but take into account the uncertainty in the q^ dependence of the axial form factor encoded in the 
axial mass Ma 

As explained above, our treatment of the reactions on nuclei relies on the local density approximation. For the 
distribution of protons in we take the empirical harmonic oscillator parametrization, and the parameter errors, 
from Ref. [s^. For the distributions of neutrons, harmonic oscillator parametrizations and relative errors as in Ref. 
are adopted but with central values from Ref. [s^. As the in-medium A(1232) broadening plays an important role, 
we have gauged its uncertainty with a 10% relative error. One should point out that our error estimate for nuclear 
effects does not account for possible multi-nucleon mechanisms that are not part of the model developed in Ref. [l^. 


^ With the helicity amplitudes -A 1/2 ^ 3/2 ^_= 0 from MAID, one obtains a branching ratio of r(A —>• A^7)/rtot(A) = 0.6%, 

consistent with the PDG estimate of 0.55-0.65% 1^ . 
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TABLE I. Error budget. 


Quantity 

Value 

Source 

Ma 

1.016 ± 0.026 GeV 

m 


1.00 ±0.11 

[31J 

Maa 

0.93 ± 0.07 GeV 


Ai/2 

(-140 ±6)10"® GeV- 1/2 

[32,^ 

A 3/2 

(-265 ±5)10-3 GeV-V2 

[32,^ 

dp 

1.833 ± 0.014 fm 

[3^ 

(y.p 

1.544 ± 0.001 fm 


dn 

1.815 ±0.014 fm 

[35,^ 

Oin 

1.529 ± 0.001 fm 

[35,^ 

(ImSA)r 

r = 1.0 ±0.1 



D. Photon events at SK 


Once the differential cross sections for the components of the detector (H 2 O) are established, it is straightforward 
to obtain the number of NC 7 events for a given photon energy and direction with respect to the neutrino beam, 


dN 


dEjd cos 


= NpOT E E Nt / dE,4>i{E,) 


dcrit{E^) 


l — t—p, 

Here, the total numbers of protons and nuclei in the SK inner detector are 


dEyd cos 


(13) 


= —MNa = -MNa , 

P 18 A g 


Nwo = —M— = —MNa ■ 
^ 18 16 18 


(14) 


where M = 2.25 x 10^° grams is the fiducial mass, and Na^ the Avogadro number. Our estimate is for the recent 
T2K Ue wpearance analysis, corresponding to a total number of protons on target (POT) A^pot = 6.57 x 10^° in 
y mode [a- The flux of the off-axis neutrino beam from Tokai has a narrow peak with median energy of 630 MeV 
at SK (see Fig. [3|). We neglect its contribution above E^ = “i GeV. In spite of the rather long tail, a sizable 
contribution of the E,^ > 3 GeV region would require a considerably large cross section at high energies, which we do 
not expect (see the discussion in Ref. i). The negative result in the single photon search performed by the NOMAD 
experiment [38l| . with an average energy of the neutrino flux of E^, ~ 25 GeV, is in line with our assumption. As the 
NG interaction is flavor independent, the composition of the beam after oscillations can be ignored. 



FIG. 3. (Color online) The T2K flux prediction at the SK detector (without oscillations). It is shown only below 5 GeV 
although it is simulated up to 30 GeV. 
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It will be instructive to consider also the (non-observable) neutrino-energy event distribution, which can be easily 
related to the integrated cross section 


dN 

dEy 


-/VpOT E E Nt(j)i{Ev)ait{Ev) ■ 

t—p, 


(15) 


III. RESULTS 

The photon energy and angular distributions of the NCy events at SK are shown in Fig. |4] for the different contri¬ 
butions, i.e. NCy on the two protons of H 2 O and on (coherent and incoherent). The contributions of the and 
Pfi components of the flux are displayed. The incoherent reaction is the largest and peaks at E^ ^ 200 — 300 MeV, 
reflecting the importance of the A(1232). The yield of nucleons and the coherent channel, both similar in size, is 
smaller but still important. The yield is quite small while other flux components are totally negligible. The 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



-1.0 -0.5 0.0 0.5 1.0 


(GeV) 


cosG.^ 


FIG. 4. (Color online) Photon energy (left) and angular (right) distributions of NC 7 events at SK from the T2K flux. The 
curves labeled as “p”, “inc” and “coh” stand for the contributions of the iz(p) —H 2 , i'(t')—incoherent and coherent reactions, 
respectively. Upper (lower) plots are for the (P^) component of the flux. 


angular distributions are forward peaked, particularly for the coherent reaction. The latter has larger incidence than 
the nucleon channel in the forward direction. On the other hand, the angular dependence for the incoherent reaction 
induced by antineutrinos is softer than the neutrino one and peaks around cos 9^ =0.7. Similar features in the energy 
and angular distributions were obtained for the MiniBooNE detector [2lj. 

Summing over all bins in the histograms above, one finds that the total number of NC 7 events is 

M = -p = (0.412 ± 0.049) -b (0.015 ± 0.002) = 0.427 ± 0.050. (16) 

The errors correspond to a standard 68 % confidence level (CL) and are dominated by the uncertainty in (7^(0). This 
is a small quantity compared to the 28 e-like events detected at SK @ 0 but can be relevant in future attempts to 
measure 5cp- This result becomes even more significant when compared to the NEUT 5.1.4.2 equivalent figure of 

AAneut = A7j^eut + A/'neut = 0-209 + 0.008 = 0.217. (17) 


^ The comparison is only indicative due to the lack of efficiency correction in our estimate. 





































7 


Indeed, using the NC 7 cross section model of Ref. [T^ we predict nearly twice more events than NEUT 5.1.4.2, on 
which the T2K estimate is based. In NEUT, NC 7 interactions proceed via baryon resonance excitation (predominantly 
A(1232)) followed by radiative decay. Resonance production is implemented according to the Rein-Sehgal original 
model but with the parameter mA = 1.21 GeV in the transition form factors of Eq. (3.12) in that reference. 

In view of this discrepancy, we have performed more detailed comparisons with the NEUT 5.1.4.2 NC 7 prediction 
and confronted photon energy and angular distributions. Figure[S]does not reveal significant shape differences between 
the two models although the one of Ref. [l^ predicts considerably more events at = 0.1 — 0.2 GeV and a faster 
decrease for ss 0.3 — 0.6 GeV. For the flux, the photon angular distribution from the model of Ref. [T^ is 
slightly flatter than the one from NEUT. The neutrino energy dependence of the events, displayed in Fig. [5] also shows 




(GeV) 



-1.0 -0.5 0.0 0.5 1.0 

cosG.^ 


FIG. 5. (Color online) E-y and cosOy distributions of NC 7 events at SK for the T2K experiment, obtained with the model of 
Ref. [ 1 ^ (solid line) and with NEUT 5.1.4.2 (triangles). The bands correspond to a 68 % CL based on the error budget of 
Sec. Ill Cl For a better shape comparison, NEUT results have been multiplied by a factor of 2. 


a good agreement in the shape although for the flux, the results by the model of Ref. [T^ do not show a peak 
around Ep^ ~ 1 GeV. In both models, most events come from neutrinos with energies between 0.4 and 1.2 GeV where 
the approach, choice of mechanisms and approximations of Ref. [l5l| are applicable. From the Ep^ distributions, it 
is apparent that the yield from is underestimated when only Ep^ < 3 GeV are kept because a sizable fraction of 
the Dp flux is left out. Nevertheless, this does not affect the comparisons, that have been consistently performed, nor 
the predictions for the total number of events, which are largely dominated by the Vp contribution. Indeed, in 

Eq. (fTBl) is smaller than the error in Af. 

The comparisons in Figs. [S] and [5] indicate that the disagreement is due to a discrepancy in the size (normalization) 
of the integrated cross sections in the two models. This is consistent with the comparison of the NG 7 integrated cross 
sections on ^^G from different models, as a function of Ep displayed in Fig. 9 of Ref. (40l| . 


IV. CONCLUSIONS 


The microscopic model of Ref. [1^ for single photon emission in (anti)neutrino NC interactions has been applied 
to predict the number of such events at the inner SK water Cherenkov detector, as well as their energy and angular 
distributions. With this model one can take into account not only the radiative decay of weakly excited A(1232) 
resonance on both nucleons and nuclei, but also smaller^ although relevant, contributions from nucleon pole terms 
and the coherent channel (details can be found in Ref. 1^). For a VpoT = 6.57 x 10^° we predict 0.427 ± 0.050 
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FIG. 6 . (Color online) Ei, distribution of NC 7 events at SK for the T2K experiment. Line styles have the same meaning as in 
Fig. [ 5 ] The results from the NEUT generator have been multiplied by a factor of 2. 


events without efficiency corrections; only 3% of it arises from the Dfj, contamination of the ly/j, beam. This small but 
irreducible background should be realistically estimated in order to increase the precision in the determination of 
oscillation parameters, particularly in Scp measurements. 

Remarkably, the prediction based on the model of Ref. [l^ is twice larger than the one obtained from the main 
T2K Monte Carlo generator NEUT (version 5.1.4.2). In a detailed inspection, we have found no significant 
differences in the shapes of the photon energy, photon angular and neutrino energy distributions in the two models. 
The large difference in normalization cannot be solely attributed to the lack of non-A production amplitudes or 
coherent photon emission in NEUT. The same is true for the different treatment of nuclear corrections. A closer 
inspection of the dominant A(1232) mediated mechanism on the nucleon indicates that the largest mismatch arises, 
not from the A(1232) production, common to weak pion production reactions, but from the strength of the A(1232) 
radiative decay, which is smaller in NEUT 5.1.4.2. It is also worth recalling that, as shown in Ref. [2lj, the number 
of NC 7 events at the MiniBooNE detector predicted by the same microscopic model used here is consistent with the 
MiniBooNE in situ estimate obtained with the NUANCE generator [4l| tuned to the NCtt^ measurement [4^. 
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